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■ Abstract 

> : 

ly-) , A general theory of the interaction of the quantized electromagnetic field 

CN ' with atoms in the presence of dispersing and absorbing dielectric bodies of 

• given Kramers-Kronig consistent permittivities is developed. It is based on 

| a source-quantity representation of the electromagnetic field, in which the 

' electromagnetic-field operators are expressed in terms of a continuous set of 

fundamental bosonic fields via the Green tensor of the classical problem. In- 
'"qJ troducing scalar and vector potentials, the formalism is extended in order to 

i ' include in the theory the interaction of the quantized electromagnetic field 

with additional atoms. Both the minimal-coupling scheme and the multipolar- 
coupling scheme are considered. The theory replaces the standard concept of 
qh mode decomposition which fails for complex permittivities. It enables us to 

treat the effects of dispersion and absorption in a consistent way and to give 
a unified approach to the atom-field interaction, without any restriction to a 
particular interaction regime in a particular frequency range. All relevant in- 
formation about the dielectric bodies such as form and intrinsic dispersion and 
absorption is contained in the Green tensor. The application of the theory to 
the spontaneous decay of an excited atom in the presence of dispersing and 
absorbing bodies is addressed. 
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1 Introduction 

Optical instruments such as beam splitters or cavities are more or less complicated 
macroscopic material bodies, whose use in experiments requires careful examination 
with regard to their action on the light under study. In quantum optics an impor- 
tant consideration is the influence of the presence of material bodies on the quantum 
statistics of the light. It is therefore necessary to take account of the presence of 
material bodies when considering the interaction of quantized light with atomic sys- 
tems. In principle, such bodies could be included as a part of the matter to which 
the radiation field is coupled and treated microscopically. However there is a class of 
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material bodies whose action can be included in the quantum theory exactly, namely 
dielectric bodies that respond linearly to the electromagnetic field. 

The quantum theory of radiation in the presence of dielectric media has been 
studied over a long period. Commonly, dielectric matter is characterized by the per- 
mittivity, which describes the response of the matter to the electric field. Quantization 
of the electromagnetic field in dielectrics with real and frequency-independent per- 
mittivity has been treated extensively 0, |, |, |, g, |, [7], |, |, 0, |TJ, g|, [T3|. In the 



same context, dispersive dielectrics have been considered |jT^, [15], [16], [FT], |Tj| [19], [20| . 
However, it is well known that the permittivity is a complex function of frequency 
which has to satisfy the Kramers-Kronig relations which state that the real part 
of the permittivity (responsible for dispersion) and the imaginary part (responsi- 
ble for absorption) are necessarily connected with each other. Hence, any quantum 
theory that is based on the assumption of a real permittivity can only be valid for 
narrow-bandwidth fields far from medium resonances where absorption can safely be 
disregarded. 

A systematic and quantum-theoretically consistent approach to the problem was 
developed by Huttner and Barnett |H|] on the basis of the microscopic Hopfield model 
of a bulk dielectric [p2^] . They performed an explicit Fano-type diagonalization J23] 
of a Hamiltonian consisting of the electromagnetic field, a (harmonic-oscillator) po- 
larization field representing the dielectric matter, and a continuous set of (harmonic- 
oscillator) reservoir variables accounting for absorption. The resulting expression for 
the vector potential could also be written in terms of the Green function of the 
classical scattering problem as was realized by Gruner and Welsch |24| , 25 1. That, 
in fact, makes it possible to perform the quantization of the electromagnetic field 
in the presence of arbitrary dielectric bodies of phenomenologically given Kramers- 
Kronig consistent permittivities, without referring to specific microscopic models of 
the bodies, which are hard to establish for general systems [^6], [27 |. 

Having quantized the electromagnetic field, the question arises of how to include 
in the theory the interaction of the medium-assisted field with atomic systems in 
order to study the influence of dielectric bodies on quantum-optical processes such 
as spontaneous decay of an excited atom near material bodies or in micro-cavities. It 
is well known that the spontaneous decay can be strongly modified when the atom 
is placed inside a high-quality micro-cavity. Recent progress in constructing micro- 
cavities has rendered it possible to approach the ultimate quality level determined 
by intrinsic material losses, so that the question of the influence of absorbing ma- 
terial has been of increasing interest. In what follows we develop a basic-theoretical 
concept for treating the interaction of the electromagnetic field with atomic systems 
in the presence of dielectric bodies of (phenomenologically) given arbitrary complex 
permittivities. 

The paper is organized as follows. In Sec. |2| the quantization of the electromagnetic 
field on a dielectric background of complex permittivity is outlined. The interaction 
of the medium-assisted electromagnetic field with atomic systems is studied in Sec. |3] 
and the minimal-coupling and multipolar-coupling Hamiltonians are derived. Finally, 
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a summary and some concluding remarks are given in Sec. 

2 Quantization of the electromagnetic field in the 
presence of dielectric media 

As already mentioned, dielectric matter plays an important role in optics, because 
(passive) optical instruments are typically composed of dielectrics. In classical electro- 
dynamics, dielectric matter is commonly described in terms of a phenomenologically 
introduced dielectric susceptibility (or permittivity). This concept has the benefit of 
being universally valid, because it uses only general physical properties, without the 
need of involved ab initio calculations. 

2.1 Classical basic equations 

The phenomenological Maxwell equations of the electromagnetic field in the presence 
of dielectric bodies but without additional charge and current densities read 

V-B(r) = 0, (1) 

V x E(r) + B(r) = 0, (2) 

V ■ D(r) = 0, (3) 

V x H(r) - D(r) = 0, (4) 

where the displacement field D is related to the electric field E and the polarization 
field P according to 

D(r) = e E(r) + P(r), (5) 
and for nonmagnetic matter it may be assumed that 

H(r) = -B(r). (6) 

Ho 

Let us restrict our attention to isotropic but arbitrarily inhomogeneous media and 
assume that the polarization linearly and locally responds to the electric field. In this 
case, the most general relation between the polarization and the electric field which 
is in agreement with the causality principle and the dissipation-fluctuation theorem 
is 

oo 

P(r,t)=e J dr X (r,r)E(r,t-r)+P n (r,t), (7) 

o 

where x( r ; T ) is the dielectric susceptibility as a function of space and time, and P n 
is the (noise) polarization associated with absorption. 
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Substitution of this expression into Eq. (|5|) together with Fourier transformation 
converts this equation to 

D(r, u) = s s(r, w)E(r, cu) + P n (r, u), (8) 

where 

oo 

e(r,o;) = l+ / dr X (r,r)e^ (9) 



is the (relative) permittivity, and the Maxwell equations (§])-(§J) read in the Fourier 
domain as 

V-B(r,w) = 0, (10) 

V x E(r, a) = iwB(r, u), (11) 
£ V-£(r,o;)E(r,a;) = pjr,u), (12) 

V x B(r, to) + i-ge(r, w)E(r, w) = //ojfr, w) (13) 

(c~ 2 = eoyUo )-Q Here we have introduced the noise charge density 

P n (r,u) = -V-P n (r,u) (14) 

and the noise current density 

3 n (r,u) = -iuP n (r,uj), (15) 
which obey the continuity equation 

V • j n (r,u})=iup D (r,cj). (16) 
According to Eq. the permittivity e(r,uj) is a complex function of frequency, 

e(r,uj) = e'(r,uj)+ie"(r,uj). (17) 

The real and imaginary parts, which are responsible for dispersion and absorption 
respectively, are uniquely related to each other through the Kramers-Kronig relations 

7T J UJ — UJ 

e (r,uj) = / dc^ 19 

7T J UJ — UJ 

(V, principal value). Further, e(r, u) as a function of complex w satisfies the relation 

e(r,-a;*)=£*(r,u;) (20) 



1 Here and in the following the Fourier transform F_(ui) of a real function F(t) is defined according 
to the relation F(t) = f™ dw F(oj)e- iut + c.c. 
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and is holomorphic in the upper complex half-plane without zeros. In particular, it 
approaches unity in the high-frequency limit, i.e, e(r, uj) — *■ 1 if \u\ — ► oo. 

The Maxwell equations ([TTJ) and flT3| ) imply that E(r, uj) obeys the partial differ- 
ential equation 

uj 2 

V x V x E(r,w) - — e(r, w)Ej>,cj) = ^/i j n (r,w), (21) 
whose solution can be represented in the form 

E(r,oj) = i/i Q uj / dVG(r,r» • j n (r», (22) 

where the Green tensor G(r, r', cu) has to be determined from the equation 

uj 2 

V x V x G(r, r', uj) - — e(r, w)G(r, r', w) = <5(r - r') (23) 

<r 

together with the boundary condition at infinity. In Cartesian coordinates, Eq.(f2~3|) 
reads 

G mi (r,r>)=<^(r-r') (24) 



(dld r m - 5 im A r ) - 5 im —e(r,uj) 

c A 

{d\ = d/dxi), where over repeated vector-component indices is summed. The Green 
tensor has the properties that 

Gij(r,r',uj)* = Gij(r,r',-uj), (25) 

G ji (i f ,T,u) = G ij (T,i f ,u), (26) 

and 2 

J d 3 s ^e"(s, uj) G lk (r, s, uj)G* k {v' , s, uj) = Im G -(r, r', uj). (27) 

The property (1251) is a direct consequence of the corresponding relation (|20|) for the 



permittivity, Eq. (|26|) is the so-called reciprocity relation, and (|27|) is proved in |27| . 

The Fourier components of the magnetic induction, B(r, uj), and the displacement 
field, D(r, uj), are directly related to the Fourier components of the electric field, 
E(r,w), 

B(r, w) = {iujY 1 V x E(r, cu), (28) 

D(r, w) = (/Wr 1 V x V x E(r, w) (29) 

[see Eqs. (PI), ©, (0), and (|2l|) 1, and E(r, uj) is determined, according to Eq. ([22]), 
by j n (r, w). The continuous set of (complex) fields j n (r, uj) [or, equivalently P n (r, uj)] 
can therefore be regarded as playing the role of the set of dynamical variables of 
the system composed of the electromagnetic field and the medium (including the 
dissipative system). For the following it is convenient to split off some factor from 
P n (r, uj) and to define the fundamental dynamical variables f(r,u) as follows: 



Pjr.,;) IX ^s"{v,uj)i{v,uj). (30) 
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2.2 Field quantization 

The transition from classical to quantum theory now consists in the replacement of 
the classical fields f(r, uj) and f*(r, uj) by the operator- valued bosonic fields f(r, uj) 
and V(r,uj) respectively, which are associated with the elementary excitations of the 
composed system in linear approximation. Thus the commutation relations are 

[f k (r,u),fl(r',u')} = 6 kk ,5(r-r')5(uj-uj'), (31) 

[/ fc (r,u;),A,(rV)] =0, (32) 
and the Hamiltonian of the composed system is 

oo 

H = J d 3 r J duhuf j (r,u)-i(r,u) (33) 
o 

Recalling Eqs. (JT^) and fl37J) and replacing E(r,w) [Eq. @], B(r,u) [Eq. 
and D(r, uj) [Eq. (p9|)] by the quantum-mechanical operators, we find that 



|r )W ) = i*/-^ /dVv^MG^r'.^-flrV), (34) 



7re c 2 

B(r, w) = (iuj)- 1 V x E(r, w), (35) 

and 

D(r, w) = (fiouj 2 )-^ x V x E(r, w), (36) 

from which the electromagnetic field operators in the Schrodinger picture are obtained 
by integration over uj: 

oo 

E(r) = E m (r) = / E(r, to) + H.c., (37) 



B(r) = / dwB(r,w) +H.c, (3* 



o 



oo 

D(r) = D m (r) = £>£(r) = / do; D(r, uj) + H.c. (39) 



o 



In this way, the electromagnetic field is expressed in terms of the classical Green tensor 
Gr(r, r', to) satisfying the generalized Helmholtz equation (|23|) and the continuum of 



2 The longitudinal (F") and transverse (F- 1 ) parts of a vector field F are defined by F^- L ^(r) = 
J d 3 r' (r — r')F(r'), with i'"^(r) and S^' (r) being the longitudinal and transverse tensor- valued 
(5-functions respectively. 
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the fundamental bosonic field variables f(r, lu) [and V(r, u)]. All the information 
about the dielectric matter (such as its formation in space and its dispersive and 
absorptive properties) is contained [via the permittivity e(r,uj)} in the Green tensor 
of the classical problem. Eqs. fl37|) - (0), together with Eqs. (|4j) - (|36|) can be 
considered as the generalization of the familiar mode decomposition. 

The quantization scheme meets the basic requirements of quantum electrodynam- 
ics. So it can be shown by using very general properties of the permittivity and the 
Green tensor that E and B satisfy the correct (equal-time) commutation relations 



27] 



[E k (r),E k ,(r')] = = [B k (r), B k ,(r')] , (40) 
[e E k (r),B k ,(r')] = -ih e km d[8{v - r'). (41) 

Obviously, the electromagnetic field operators in the Heisenberg picture satisfy the 
Maxwell equations (|l]) - with the time derivative of any operator Q being given 
by 

$ = (ih)- 1 [Q,H], (42) 

where H is the Hamiltonian (|33|), 

Further, scalar (jp) and vector (A) potentials can be introduced and expressed in 
terms of the fundamental bosonic fields. In particular, the potentials in the Coulomb 
gauge are defined by 

-Vy3(r) = E"(r), (43) 

oo 

A(r) = J d^A(r,w) + H.c. , (44) 
o 

where 

A(r, u) = (icu)" 1 E ± (r, u) (45) 
The canonically conjugated momentum field with respect to A(r) is 

oo 

n(r) = -is J du wA(r, to) + H.c, (46) 
o 

and it is not difficult to verify that n = — ^oE^, V x A = B, and —A — Vy3(r) = E. 
In addition, A and II satisfy the well-known commutation relations 

[A k (r),A k ,(r')] = = [n fc (r),n fc , (!■')], (47) 
[A k (r),il k ,(r>)]=ihSi k ,(v-v'). (48) 
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3 Interaction of the medium-assisted quantized 
electromagnetic field with charged particles 

The interaction of the quantized electromagnetic field with atoms placed inside a 
dielectric medium or near dielectric bodies can be strongly influenced by the dielectric 
medium. A well-known example is the dependence of the spontaneous decay rate of 
an excited atom on the properties of an dielectric environment. In order to study 
such and related phenomena, the Hamiltonian ( ]3"B| ) must be supplemented with the 
the Hamiltonian of additional charged particles and their interaction energy with the 
medium-assisted electromagnetic field. 

3.1 The minimal-coupling Hamiltonian 

Applying the minimal- coupling scheme, we may write the the total Hamiltonian in 
the form 

oo 

H = / d 3 r / duhw&(r,uj)-i(T,uj) + V- [p a - g a A(r a ) 

J J ^ 2m Q L 

o 

+| j d 3 rp A (r)^ A (r) + J d 3 r p A {r)<p{v), (49) 

where r a is the position operator and p Q is the canonical momentum operator of 
the ath (non-relativistic) particle of charge q a and mass m a . The Hamiltonian ( f49"|) 
consists of four terms. The first term is the energy of the electromagnetic field and 
the medium (including the dissipative system), as introduced in Eq. (p3|). The second 
term is the kinetic energy of the charged particles, and the third term is their Coulomb 
energy, where the corresponding scalar potential 0\ is given by 

fPA(r) = / dV ^ , (50) 
J 47re |r-r'| 

with 

PA(r) = ^2qJ(r-r a ) (51) 

a 

being the charge density. The last term is the Coulomb energy of interaction of the 
particles with the medium. From Eq. fl4~9|) it follows that the interaction Hamiltonian 

reads 

H-mt = P Q ~ |<?aA(r Q ) g Q A(r a ) + / d 3 r p A (r)0(r). (52) 

Tfl a I J J 

a 

Note that in Eq. ([52]) the scalar potential <p and the vector potential A must be 
thought of as being expressed, on using Eqs. ([43]), (p|) , (^5|) together with Eqs. fl34|) 
and (P7|), in terms of the fundamental fields f(r, u) and f'(r, ui). 
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In a straightforward but somewhat lengthy calculation it can be shown that both 
the operator-valued Maxwell equations 



V -B(r) = 0, 

V x E(r) + B(r) = 0, 
V-D(r) = p A (r), 

VxH(r)-D(r)=j A (r), 
and the operator-valued Newtonian equation of motion 



1 

m a 
i i 



p a - q a A(r c 



2 , r a x B(r a ) - B(r a ) x r Q 



are fulfilled. In Eq. (|56|), the atomic current density j a (r) reads 



(53) 

(54) 
(55) 

(56) 

(57) 
(58) 

(59) 



where [ , ]+ denotes the anticommutator. Note that compared with Eqs. ( p7|) and 
(|39"|), the electric and displacement fields now contain additional longitudinal parts 
that result from the charge distribution />a(i*), i.e., 



E(r) = E m (r) - V^ A (r) 



D(r) = D m (r) - e o V0 A (r) 



du;E(r,u;) + H.c. 



dtuD(r,cj) + H.c. 



- V^ A (r) , 



(60) 



-e o V0 A (r). (61) 



The Maxwell equations (|53|) and (|55| ) simply result from the definition of the field 
operators B(r) [Eqs. (^j and (|§] and D(r) [Eqs. (||), (0), and (§J)] respectively. 
The other Maxwell equations ([54] ) and (|56|) and the Newtonian equation of motion 
( p?|) and (|58D follow from the Heisenberg equation of motion, Eq. (2.37), with the 
Hamiltonian H from Eq. (|49|). 



3.2 The multipolar-coupling Hamiltonian 

In the minimal-coupling scheme, the interaction Hamiltonian (|52|) is expressed in 
terms of the potentials of the medium-assisted electromagnetic field. With regard to 
(localized) atomic systems (atoms, molecules etc.) the interaction is commonly de- 
sired to be treated in terms of the electromagnetic field strengths and the atomic 
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polarization and magnetization. This can be achieved by means of a unitary trans- 
formation. 

Let us consider an atomic system localized at position r A and introduce the atomic 
polarization 

i 



P A (r) = Qa ( r Q - r A) / dA 5[r - r A - A (r Q - r A )] 

o 

so that the charge density (|5l|) can be rewritten as 

/5 A (r) = ^g a <5(r-r A ) - V-P A (r). 



(62) 



(63) 



In order to perform the transition from the minimal-coupling scheme to the 
multipolar-coupling scheme, we apply to the variables the unitary operator 



U = exp 



- / d 3 rP A (r)-A(r) 



Pa = Up a ll j = p Q - g a A(r a ) - J d 3 rn Q (r) x B(r), 



f(r,u)-L —e"(r,u)^ / d : V Pitr') • G' (r. r. ^) 

h V 7T£o c 



where the abbreviation 



n Q (r) = q a (r Q -r A ) / dA A 8[r - r A - A (r Q -r A )] 



(64) 



It is not difficult to prove that the following transformation rules are valid: 

r' a = m a U ] = r a , (65) 



(66) 



(67) 



(68) 



has been used. Employing equations ( |65"D - ( |67|) ) we can express the Hamiltonian H 
in Eq. as the unitary transform of a new Hamiltonian 7i, 

H = UHU\ (69) 

where 

oo 

H = J d 3 r y &u%ujf\r,u)-f(r,uj) + p a + Jd 3 rh a (r) x B(r 



+ l - J d 3 rp A (r)</$ A (r) + ±- I d 3 rP|(r) • P|(r) 



d 3 rP^(r) • E m (r) + / d 3 r p A (r)£(r). 



(70) 
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In particular when the charged particles form a neutral atomic system (%2 a q Q = 0), 
then Eq. ([?(]) takes the form of 



oo 

U = J d 3 r J dcj/^f t (r,cu)T(r,cu) + ^^- p a + Jd 3 rh a i 



1 2 



r) x Bfr) 



1 

27 



f— / d 3 rP A (r)-P A (r) 



d 3 rP A (r)-E m (r). 



(71) 



Note that the last term on the right-hand side in Eq. ( [fl]) describes the interac- 
tion between the atomic polarization and the whole medium-assisted electric field. 
Introducing the polarization of the medium 



D — eqE — D m — e E m — — e E rj 



(72) 



[see Eqs. ( |60D and (|6T|| and recalling that = D -1 , we may rewrite the last term in 
the transformed Hamiltonian ([71]) to obtain 



- /d 3 rP A (r)-E m (r) 



1 



to 



d 3 rP A (r) ■ D ± (r) + - / d 3 rP A (r) • P m (r). (73) 



Combining Eqs. (p5|), (|7TD, and Q73|), we express in the original Hamiltonian the old 
variables in terms of the new ones. The result is the multipolar-coupling Hamiltonian 



oo 

H = J d 3 r J duhuf'\r,uj)-{'{r,uj) + J2^- Pa + y"d 3 rn 



1 2 



It) x Bfr) 



1 

2^ 



d 3 rP A (r) • P A (r) 



to 



d 3 rP A (r) •D'- L (r) + 



d 3 rP A (r)-P m (r) 



where 



D- L (r)=D J -(r)+Pi(r). 
From Eq. ([74]) the interaction Hamiltonian is seen to be 

Hint' = -- / d 3 rP A (r) • D'^r) + - [ d 3 rP A (r) • P m (r 



(74) 
(75) 



+ Pl>|d 3 rn a (r)xB(r) + ^ ^— J d 3 rn Q (r) x B(r) 



(76) 



The first term on the right-hand side in Eq. (|7E| ) describes the interaction of the 
polarization of the atomic system with the transverse part of the overall displacement 
field [cf. Eq. ([75|) 1 . The second term is a contact term between the medium polarization 
and the polarization of the atomic system. The last two terms refer to magnetic 
interactions. 
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4 Summary and outlook 

We have developed a general theory of the interaction of the quantized electromag- 
netic field with atoms in the presence of dispersing and absorbing dielectric bodies 
of given Kramers-Kronig consistent permittivities. The concept is based on a source- 
quantity representation of the electromagnetic field, in which the electromagnetic- 
field operators are expressed in terms of a continuous set of fundamental bosonic 
fields via the Green tensor of the classical problem. The theory, which is a natural 
extension of the standard concept of mode decomposition, gives a unified approach 
to the atom-field interaction, without any restriction to a particular frequency range. 

The formalism has been applied quite recently to the problem of the spontaneous 
decay of an excited atom in the presence of dielectric bodies f28 |. In particular, it 



has been shown that the temporal evolution of the atomic upper-state-probability 
amplitude C u {t) obeys the integral equation 



t 



C u (t) = 1 + J dt'K(t-t')C u (t') } (77) 
o 

where the kernel function is determined by the Green tensor (at the atomic position 
ta) as follows: 

oo 

K(t - t') = [ dc Im ^ (rA ' r t U;) \e<™ *W - ll (78) 

o 

(wa, transition frequency; //j, transition dipole moment). It is worth noting that the 
integral equation ( |77| ) applies to the spontaneous decay of an atom in the presence of 
an arbitrary configuration of dispersing and absorbing dielectric bodies. All the matter 
parameters that are relevant for the atomic evolution are contained, via the Green 
tensor, in the kernel function (f78| ). It should be pointed out that the Green tensor has 
been available for a large variety of configurations such as planarly, spherically, and 
cylindrically multilayered media |29[. A first evaluation of the integral equation ([77]) 
for the case of the atom being placed at the centre of a spherical micro-cavity whose 
wall is modeled by a band-gap dielectric of Lorentz type has been given in reference 



28 to which the reader is referred for details. 
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